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ANALYSIS  OF  TURBULENT  HEAT  TRANSFER,  MASS  TRANSFER,  AND 
FRICTION  IN  SMDOTH  TUBES  AT  HIGH  PRANDTL 
AND  SCHMIDT  NUMBERS 
By  Robert  G.  Deissler 


SUMMARY 

The  expression  for  eddy  diffusivity  from  a previous  analysis  was 
modified  in  order  to  account  for  the  effect  of  kinematic  viscosity  in 
reducing  the  turbulence  in  the  region  close  to  a wall.  By  using  the 
modified  expression,  good  agreement  was  obtained  between  predicted  and 
experimental  results  for  heat  and  mass  transfer  at  Prandtl  and  Schmidt 
numbers  between  0.5  and  3000.  The  effects  of  length-to-diameter  ratio 
and  of  variable  viscosity  were  also  investigated  for  a wide  range  of 
Prandtl  numbers. 


INTRODUCTION 

Most  of  the  existing  analyses  for  turbulent  heat  and  mass  transfer 
are  adequate  only  for  Prandtl  and  Schmidt  numbers  on  the  order  of  1 or 
less.  For  instance,  the  analysis  given  in  reference  1,  although  adequate 
for  gases,  gives  heat-  and  mass-transfer  coefficients  for  liquids  at  high 
Prandtl  or  Schmidt  numbers  that  are  higher  than  those  obtained  experimen- 
tally. The  difference  between  the  experimental  values  and  the  values  ob- 
tained by  the  method  in  reference  1 increases  as  the  Prandtl  or  Schmidt 
number  increases.  Coefficients  obtained  from  the  von  Karman  analysis 
(ref.  2)  at  high  Prandtl  or  Schmidt  numbers  are  lower  than  the  experi- 
mental values.  Rannie's  analysis  (ref.  3)  gives  coefficients  which  are 
in  somewhat  better  agreement  with  the  data  than  either  of  these  analyses, 
but  the  coefficients  are  again  inaccurate  at  very  high  prandtl  or  Schmidt 
numbers.  The  analysis  of  reference  4 agrees  with  data  at  Prandtl  or 
Schmidt  numbers  of  1 and  at  very  high  Prandtl  or  Schmidt  numbers,  but  the 
coefficients  are  somewhat  low  at  intermediate  values  of  these  numbers. 

In  reference  5,  which  was  published  since  the  present  investigation  was 
initiated,  good  agreement  was  obtained  with  mass -transfer  data  for 
Prandtl  and  Schmidt  numbers  between  0.5  and  3000  by  introducing  an 
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appropriate  amount  of  turbulence  into  the  laminar  sublayer.  In  all 
these  analyses,  except  those  in  references  1 and  3,  the  properties  were 
constant  and  the  flow  fully  developed.  The  relations  among  heat  trans- 
fer, mass  transfer,  and  fluid  friction  are  discussed  in  reference  6. 

The  inadequacy  of  most  of  the  previous  analyses  at  high  Prandtl  and 
Schmidt  numbers  is  principally  caused  by  the  expressions  used  for  the 
eddy  diffusivity  in  the  region  very  close  to  the  wall.  This  region  is 
important  because  of  the  extremely  large  temperature  or  concentration 
gradients  in  that  region  at  high  Prandtl  or  Schmidt  numbers  (ref.  7). 

In  the  analysis  given  herein,  which  was  made  at  the  NACA  Lewis  labora- 
tory, the  expression  for  eddy  diffusivity  given  in  reference  1 is 
modified  in  order  to  account  for  the  effect  of  kinematic  viscosity  in 
reducing  the  turbulence  in  the  region  close  to  the  wall.  The  effects 
of  variable  viscosity  and  of  length-to -diameter  ratio  are  also 
investigated. 


BASIC  EQUATIONS 

'For  obtaining  the  velocity,  temperature,  and  concentration  distri- 
butions in  a tube  with  turbulent  flow,  the  differential  equations  for 
shear  stress,  heat  transfer,  and  mass  transfer  can  be  written  as  fol- 
lows (symbols  are  defined  in  the  appendix): 


_ du  , du 

T * 11  Sy  + pt  3y 

(1) 

, dt  dt 

= "k  dy  ‘ pscPeb  dy 

(2) 

m _ -v  dC  , dC 
dy  h dy 

(3) 

where  the  values  for  e and  are  dependent  on  the  amount  and  kind 

of  turbulent  mixing  at  a point.  The  eddy  diffusivities  for  heat  and 
mass  transfer  are  equal  inasmuch  as  both  processes  are  governed  by  the 
same  differential  equation  if  aerodynamic  heating  is  neglected  (diffusion 
equation).  On  the  other  hand,  the  ratio  e^A  = a must  be  determined 

experimentally  or  theoretically  inasmuch  as  the  equation  of  motion  for 
a fluid  contains  terms  which  are  not  present  in  the  diffusion  equation. 
Equations  (l)  to  (3)  can  be 


written  in  dimensionless  form  as 


(_P_  + P_  c \ du+ 
VQ  PO  Vq/p0  J dy+ 


l 

\ 


\ 
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q. 

‘Jo 


/ k 1 p cp  e \ dt+ 

1*0  Pr0  Pq  cp>0  Pq/pq  ) dy+ 


i - (±  1 + a -4_\ 

“0  \\)  SCq  Mq/Po  / dy+ 


(5) 


(6) 


The  variation  of  properties  in  these  equations  might  be  caused  by 
either  radial  variation  of  temperature  or  of  concentration  of  the 
diffusing  substance. 


EXPRESSIONS  FOR  EDDY  DIFFUSIVITY 


In  order  to  make  practical  use  of  equations  (4),  (5),  and  (6),  the 
eddy  diffusivity  e must  be  evaluated  for  each  portion  of  the  flow. 
After  consideration  of  the  various  factors  on  which  e might  depend, 
the  following  functional  relation  is  assumed: 


e m 


e 


u du 
PJ  dy' 


d2u  d3u 
dy2*  dy3’ 


(7) 


(The  quantities  u and  p must  occur  together  inasmuch  as  they  are 
the  only  quantities  in  eq.  (7)  containing  mass  as  a dimension.)  It  is 
assumed  in  reference  8 that  in  the  region  at  a distance  from  the  wall 

e = e(du/dy,  d2u/dy2)  (Karman's  as sumpt ion ) , and  in  the  region  close  to 
the  wall  e ■ e(u,y)  - n^uy.  In  both  of  these  expressions  the  possible 
effect  of  kinematic  viscosity  u/p  is  neglected.  It  appears  from  heat- 
and  mass-transfer  data  at  high  Prandtl  and  Schmidt  numbers,  however, 
that  the  effect  of  p/p  cannot  be  neglected  in  the  region  very  close  to 

the  wall  (y+ < 5).  When  the  previous  expressions  are  used  to  compute 
heat  and  mass  transfer  at  high  Prandtl  and  Schmidt  numbers,  the  coeffi- 
cients are  considerably  too  high  compared  with  the  experimental  values. 
It  might  be  expected  that,  in  the  region  very  close  to  the  wall  where 
the  turbulence  level  is  low,  the  effect  of  kinematic  viscosity  would  be 
important  inasmuch  as  the  viscous  effects  would  be  of  the  same  order  of 
magnitude  as  the  inertia  effects.  Therefore,  for  the  region  close  to 
the  wall,  e is  written  as 


\ 


\ 


e = e(u,y,u/p) 


(8) 
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As  in  reference  8,  the  effect  of  the  derivatives  is  neglected  close  to 
the  wall  because  the  flow  is  very  nearly  laminar  in  that  region.  The 
first  derivative  approaches  the  value  u/y  and  hence  may  be  omitted 
since  u and  y already  appear  in  the  functional  relation.  The  second 
derivative  becomes  very  nearly  zero  as  the  wall  is  approached. 

From  dimensional  analysis,  equation  (8)  becomes 


e 


n2uyF 


0) 


The  function  F[e'/(u/p)j  should  approach  1 as  t or  t'  increases, 
because  the  effect  of  kinematic  viscosity  becomes  negligible  at  high 
turbulence  levels.  Inasmuch  as  dimensional  analysis  cannot  determine 
the  form  of  the  function  F,  additional  assumptions  must  be  made. 

The  simplest  assumption  that  might  be  made  for  F[e'/(u/p)]  is 
that  it  equals  e'/(n/p).  This  assumption  could  be  written  in  differen- 
tial form  as 


dF  = 


(10) 


However,  equation  (10)  could  not  be  expected  to  hold  as  F approaches 
1,  because  the  change  in  F for  a given  change  in  e'/(p/p)  should 
approach  0 as  F approaches  1 (F  could  never  be  greater  than  l). 

The  simplest  multiplicative  factor  which  gives  equation  (10)  this  char- 
acteristic is  (l  - F).  The  assumption  made  for  the  variation  of  F, 
the  adequacy  of  which  will  be  checked  by  experiment,  therefore 
becomes 


dF  = d[e'/(u/p)](l  - F) 


(11) 


Separating  variables  and  integrating  from  the  wall  to  a point  in  the 
fluid  give 


lt"/(u/p) 


(12) 
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where  the  lower  limits  are  taken  as  zero,  because  the  turbulence  goes 
to  zero  at  the  wall  and  the  effect  of  kinematic  viscosity  consequently 
becomes  extremely  large  at  the  wall  (F  -*■  0).  Integration  of  equation 
(12)  results  in 

" (m/p) 

F * 1 - e (13) 


This  expression  approaches  1 for  large  values  of  e'/(M/p)*  From 
equations  (9)  and  (13 ) there  results 


e ■ 


(14) 


Equation  (14)  gives  e as  a function  of  u and  y for  the  region 
close  to  the  wall.  The  constant  n is  to  be  determined  experimentally. 


For  the  region  at  a distance  from  the  wall,  F is  usually  close  to 
1,  because  the  effect  of  kinematic  viscosity  is  small.  The  Karman 
expression  for  t,  which  neglects  the  variation  of  F,  can  usually  be 
used  in  that  region,  or 


c = x2  Jdu/dy)3 

(d2u/dy2)2 


(15) 


If  It  is  desired  to  take  the  variation  of  F into  account  in  the  region 
at  a distance  from  the  wall,  equation  (15)  becomes 


c = x2 


(du/dy)5 

(d2u/dy2)2 


. x2(du/dyf/(d2u/dy2)2 

M/P 


(16) 


No  attest  is  made  in  the  present  analysis  to  specify  the  mechanism 
by  which  the  kinematic  viscosity  reduces  the  eddy  diffusivity  in  the 
region  close  to  the  wall,  because  the  exact  mechanism  is  unknown.  Pos- 
sible mechanisms  are  these:  First,  the  kinematic  viscosity  might  act  as 

a damping  factor  to  reduce  the  turbulence  level  close  to  the  wall. 
Second,  it  might  help  to  orient  the  eddies  close  to  the  wall  by  damp- 
ing out  those  moving  perpendicular  to  the  wall  and  thus  reduce  tne 
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effective  turbulent  transfer.  Third,  it  might  act  to  produce  a partial 
turbulence}  that  is,  it  might  cause  the  flow  at  a point  to  be  laminar  a 
fraction  of  the  time.  The  actual  effect  of  kinematic  viscosity  might 
be  due  to  any,  or  all,  of  these  mechanisms. 


ANALYSIS  FOR  CONSTANT  FLUID  PROPERTIES 

In  order  to  solve  equations  (4),  (5),  and  (6),  the  following 
assumptions  are  made  in  addition  to  the  assumptions  concerning  the 
expressions  for  eddy  diffusivity  (eqs.  (14),  (15),  and  (16)): 

(1)  The  eddy  diffusivities  for  momentum  e and  heat  or  mass  trans- 
fer are  equal,  or  a = 1.  Previous  analyses  for  flow  of  gases  in 

tubes  based  on  this  assumption  yielded  heat-transfer  coefficients  that 
agree  with  experiment  (ref.  1).  At  low  Peclet  numbers  (Pe  = RePr),  a 
appears  to  be  a function  of  Peclet  number  (ref.  9)  but  is  approximately 

1 at  high  Peclet  numbers.  In  general,  the  Peclet  numbers  are  high  at  high 
Prandtl  numbers  for  turbulent  flow. 

(2)  The  variations  across  the  tube  or  boundary  layers  of  shear 
stress  x , the  heat  transfer  per  unit  area  q,  and  the  mass  transfer  per 
unit  area  m have  a negligible  effect  on  the  velocity,  temperature,  and 
concentration  distributions.  It  is  shown  in  figure  12  of  reference  1 
that  the  assumption  of  a linear  variation  of  shear  stress  and  heat  trans- 
fer across  the  boundary  layers  gives  very  nearly  the  same  velocity  and 
temperature  profiles  for  gases  as  those  obtained  by  assuming  uniform 
shear  stress  and  heat  transfer  across  the  boundary  layers  for  values 

of  6j1+  between  500  and  5000.  (The  boundary  layers  fill  the  tube  for 
fully  developed  flow  or  fully  developed  heat  or  mass  transfer.)  For 
small  values  of  &h+>  such  as  occur  very  near  the  entrance,  the  effect 
of  variable  heat  transfer  (or  mass  transfer)  is  checked  in  figure  9 of 
reference  10  and  found  to  be  negligible.  Although  these  checks  were 
made  for  gases  (Prandtl  or  Schmidt  numbers  close  to  l),  the  effect  of 
variation  of  heat  or  mass  transfer  per  unit  area  at  high  Prandtl  or 
Schmidt  numbers  would  be  even  less  because  the  temperature  or  concen- 
tration profile  becomes  flatter  as  the  Prandtl  or  Schmidt  number  is 
increased. 

(3)  The  molecular  shear-stress,  heat-transfer,  and  mass-transfer 
terms  in  the  equations  can  be  neglected  in  the  region  at  a distance  from 
the  wall  (ref.  1,  fig.  14). 

(4)  In  the  case  of  mass  transfer,  the  concentration  of  the  dif- 
fusing substance  is  small  enough  that  the  mass  transfer  does  not  appreci- 
ably change  the  velocity.  This  condition  is  generally  obtained  in  the 
case  of  evaporation  from  a wetted  wall  or  in  the  solution  of  the  wall 
material  in  a liquid. 
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Velocity,  temperature,  and  concentration  distributions.  - The  ex- 
presslon  for  c close  to  the  wall  (eq.  (14) ) can  be  written  In  dimen- 
sionless form  for  constant  properties  as 


n2u+y+ 


e 


n2u+y+> 


fl7) 


Equations  (4),  (5),  and  (6)  can  be  written  in  integral  form  for  the 
region  close  to  the  wall  with  the  preceding  assumptions  and  constant 
fluid  properties  as 


(18) 


(19) 


(20) 


For  the  region  at  a distance  from  the  wall,  equation  (15)  is  substituted 
in  equations  (4),  (5),  and  (6).  By  use  of  assumptions  (2)  and  (3), 
equation  (4)  becomes,  for  constant  fluid  properties, 


u 


L+  - 


(21) 


which  is  the  well-known  logarithmic  equation.  Division  of  equation  (4) 
by  equation  (5)  or  (6)  gives,  with  assumptions  (l),  (2),  and  (3), 


I 

i 


8 


NACA  TN  3145 


u+  - Ul+  = t+  - tx+  - C+  - (22) 

where  the  equations  are  integrated  from  y^+,  the  lowest  value  of  y+ 
for  which  the  equations  for  flow  at  a distance  from  the  wall  apply,  to 


The  variation  of  u+  with  y+  for  fully  developed  adiabatic  tur- 
bulent flow  from  references  8 and  11  is  shown  on  semi logarithmic  co- 
ordinates in  figure  1.  The  curves  corresponding  to  equations  (18)  and 
(21 ) are  also  shown  in  the  figure.  Equation  (18)  was  solved  by  numer- 
ical iteration  inasmuch  as  u+  occurs  on  both  sides  of  the  equation. 

The  values  of  the  constants  in  the  equations  are  n = 0.124,  x * 0.36, 

and  y^+  = 26,  as  determined  from  the  experimental  data.  These  values 

will  apply  also  to  the  temperature  distributions.  The  value  of  n in 
equation  (18)  differs  from  that  given  in  reference  8 because  equation 

(18)  includes  the  effect  of  kinematic  viscosity  on  e.  Figure  1 indi- 
cates good  agreement  of  equation  (18)  with  the  data  for  y+  < 26  and  of 
equation  (2l)  with  the  data  for  y+  > 26.  From  the  velocity-distribution 

data  it  is  difficult  to  tell  whether  e = n^uy  from  reference  8 or 
equation  (14)  for  e should  be  used  in  the  region  close  to  the  wall, 
inasmuch  as  both  expressions  give  results  which  agree  closely  with  the 
data.  The  tenperature  or  concentration  profiles  at  high  Prandtl  or 
Schmidt  numbers  are,  however,  much  more  sensitive  to  the  values  of  e 
very  close  to  the  wall  because  of  the  veiy  large  temperature  or  concen- 
tration gradients  in  that  region,  as  can  be  seen  in  figure  2. 

Generalized  temperature  or  concentration  distributions  calculated 
from  equations  (19),  (20),  and  (22)  are  presented  in  figure  2 on  log- 
log  coordinates.  Each  curve  represents  either  t+  or  C+  at  a given 
Prandtl  or  Schmidt  number,  as  can  be  seen  by  comparison  of  equations 

(19)  and  (20).  The  curves  indicate  that  the  temperature  or  concentration 
distributions  become  flatter  over  most  of  the  tube  radius  as  the  Prandtl 
or  Schmidt  number  increases.  From  equations  (19)  and  (20),  dt+/dy+  = Pr 
and  dC^/dy1"  = Sc  at  or  very  near  the  wall  so  that  the  slopes  of  the 
curves  at  the  wall  increase  with  Prandtl  or  Schmidt  number.  The  slopes 
of  the  curves  in  figure  2 near  the  wall  appear  equal  because  the  curves 
are  plotted  on  log-log  coordinates  (d(log  t+)/d(log  y+)  * 1 at  the  wall). 
Included  for  comparison  is  the  temperature  distribution  for  a Prandtl 
number  of  300  calculated  by  the  method  in  reference  1,  which  neglects 

the  effect  of  kinematic  viscosity  (e  = n2uy). 

The  sensitivity  of  the  temperature  or  concentration  distribution  at 
high  Prandtl  or  Schmidt  numbers  to  various  assunptions  for  the  turbulent 
transfer  in  the  region  close  to  the  wall  compared  with  that  of  the 
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velocity  distribution  indicates  that  the  region  very  close  to  the  wall 
could  be  studied  advantageously  by  measuring  temperatures  or  concen- 
trations at  high  Prandtl  or  Schmidt  numbers  rather  than  by  measuring 
velocities  in  that  region.  Some  work  along  these  lines  has  been  reported 
in  reference  5,  in  which  concentration  profiles  at  high  Schmidt  numbers 
were  measured  with  an  interferometer.  No  evidence  of  a purely  laminar 
layer  (linear  concentration  profile)  was  found  for  values  of  y+  as 
low  as  1.  This  result  is  in  agreement  with  the  assumption  in  the  pre- 
sent analysis,  where  the  turbulence  is  assumed  to  be  zero  only  at  the 
wall. 


i 

8 


Relations  among  Nusselt,  Reynolds,  and  Prandtl  or  Schmidt  numbers 
for  constant  properties.  - It  can  be  shown  from  the  definitions  of  the 
quantities  involved  that  the  Nusselt  numbers  for  heat  and  mass  transfer 
and  the  Reynolds  number  are  given  by 


Nu  S! 


2r0+Pr 

tb+ 


(23) 


Nu' 


2r0+  Sc 


(24) 


Re  = 2ub+r0+ 


(25) 


where 


(26) 


(27) 
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and 


The  Nusselt  numbers  in  these  equations  are  based  on  the  difference 
between  wall  and  bulk  temperatures  or  concentrations.  The  relation 
among  Nusselt,  Reynolds,  and  Prandtl  numbers  can  be  obtained  from  these 
equations  and  the  generalized  distributions  given  in  figures  1 and  2 
for  various  values  of  the  parameter  Tq+. 

Predicted  Nusselt  numbers  for  fully  developed  heat  or  mass  transfer 
are  plotted  against  Reynolds  number  for  various  values  of  Prandtl  or 
Schmidt  number  in  figure  3.  The  curve  for  a Prandtl  number  of  0.73 
agrees  very  closely  with  that  given  in  reference  1,  which  was  obtained 

by  using  c = n^uy  in  the  region  close  to  the  wall.  The  curves  for 
supercritical  water  at  higher  Prandtl  numbers  (Pr  » 6)  given  in  refer- 
ences 12  (fig.  7)  and  13  (fig.  18),  however,  are  higher  than  those  in 
figure  3;  therefore  the  curves  in  these  references  should  be  replaced 
by  figure  3,  although  the  values  of  the  reference  temperatures  for 
evaluating  the  fluid  properties  should  not  be  affected. 

Examination  of  the  curves  in  figure  3 indicates  that  the  slopes  of 
the  various  curves  are  approximately  equal  on  a log-log  plot.  (The 
slopes  would  be  more  nearly  equal  if  the  Peclet  number  effect  from 
ref.  9 were  included.)  This  result  justifies  the  usual  practice  in 
heat-transfer  investigations  of  writing  Nu  = f(Re,Pr)  as  f(Re)x  f(Pr) 
(usually  as  ReaPrb).  The  same  result  does  not  hold  for  very  low 
Prandtl  numbers  where  the  slopes  change  considerably. 

Comparison  of  analysis  and  experiment  for  fully  developed  heat  and 
mass  transfer.  - A comparison  between  predicted  and  experimental  results 
for  fully  developed  heat  and  mass  transfer  is  presented  in  figure  4,  in 
which  Stanton  number  is  plotted  against  Prandtl  or  Schmidt  number  for 
various  Reynolds  numbers.  The  predicted  Stanton  numbers  were  obtained 
from  figure  3 and  the  relation  St  = Nu/RePr.  The  symbols  represent 
mean  lines  through  data  for  heat  transfer  in  gases  (ref.  l)  and  in 
liquids  (refs.  14  to  19)  and  mass  transfer  by  evaporation  from  wetted 
walls  (refs.  20  to  22),  by  solution  of  the  wall  material  in  a liquid  (refs. 
23  and  24),  and  by  diffusion-controlled  electrodes  (ref.  25).  The  pre- 
dicted and  measured  values  are  in  good  agreement  over  the  entire  range 
of  Prandtl  and  Schmidt  numbers  shown  (0.5  to  3000).  The  agreement  for  a 
Reynolds  number  of  10,000  in  the  low  Prandtl  or  Schmidt  number  range 
would  be  improved  by  applying  the  Peclet  number  correction  from 
reference  9. 
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Simplified  equation  for  fully  developed  case  for  very  high  Prandtl 
or  Schmidt  numbers.  - At  very  high  Prandtl  numbers  the  essential  tern- 
perature  changes  take  place  in  the  region  very  close  to  the  wall  where 

u+  is  very  nearly  equal  to  y+.  Setting  u+  = y+  in  equation  (19), 
expanding  the  exponential  function  in  a series,  and  retaining  only  the 
first  two  terms  of  the  series  result  in 


dy+ 


jsj:  + n4(y+)4 


M) 

Integration  of  equation  (29)  and  evaluation  of  the  result  for 
give 


(29) 


-S—  rr3/* 
2/^2n 


(30) 


The  value  of  t+  at  y+  = • is  essentially  t-^+,  because  t+  is  very 

nearly  constant  except  in  the  region  very  close  to  the  wall  at  high 
Prandtl  numbers.  From  equations  (23)  and  (30 ), 


Nu  = nro+  Pr1/4  (31 ) 

or  the  Stanton  number,  in  terms  of  the  Prandtl  number  and  friction 
factor,  is 


where 


St 


2n  'Sf ~ 
a pr3/4 


Similarly,  for  mass  transfer, 


St 


2n  f\fT 

* ScV4 


(32) 


(33) 


where  n has  the  value  0.124  as  determined  in  figure  1.  The  relation 
between  Re  and  f is  given  in  figure  8 of  reference  1 or  the  curve 


\ 

\ 
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in  figure  10  in  the  present  report  for  adiabatic  flow  (5  ■ 0).  This 
relation  is,  of  course,  independent  of  Prandtl  number  for  constant 
properties.  Equation  (32)  or  (33)  is  indicated  by  the  dotted  line  in 
figure  4 and  is  seen  to  be  in  good  agreement  with  the  predicted  line 
obtained  previously  for  Pr>200. 

Comparison  of  various  analyses.  - A comparison  of  various  analyses 
is  given  in  figure  5.  It  can  be  seen  that  all  the  analyses  more  or  less 
converge  at  the  lower  Prandtl  or  Schmidt  numbers.  At  the  high  Prandtl 
or  Schmidt  numbers,  the  present  analysis  and  the  analyses  from  refer- 
ences 4 and  5 are  in  fair  agreement,  whereas  those  from  references  2 
and  3 diverge  considerably.  The  present  analysis  and  the  analysis  from 
reference  5 represent  the  experimental  data  about  equally  well. 


Heat  or  mass  transfer  in  entrance  region  for  uniform  wall  heat  or 
mass  flux,  uniform  initial  temperature  or  concentration  distribution, 
fully  developed  velocity  distribution,  and  constant  properties.  - For 
calculating  heat  or  mass  transfer  in  the  entrance  region,  it  is  assumed, 
as  in  reference  10,  that  the  effects  of  heat  or  mass  transfer  are  con- 
fined to  fluid  layers  close  to  the  surface  (boundary  layers  for  heat  or 
mass  transfer).  The  temperature  or  concentration  distributions  outside 
the  boundary  layers  are  assumed  uniform,  and  the  temperature  or  concen- 
tration is  constant  along  the  length  of  the  tube  for  the  region  outside 
the  boundary  layer.  Inside  the  boundary  layer  the  temperature  or  con- 
centration distribution  is  obtained  from  figure  2.  Integral  heat-  or 
mass-transfer  equations  are  used  for  calculating  the  thickness  of  the 
boundary  layers  for  various  distances  from  the  entrance.  It  is  shown  in 
reference  10  that  for  heat  transfer  the  relation  between  &h+  and  x/D 

for  constant  properties  is  given  by 


X 

D 


y+)<3y+ 


(34) 


A similar  equation  can  be  written  for  mass  transfer  if  the  concentration 
of  the  diffusing  material  is  small  (assumption  4): 


X 

5 = 


(35) 
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In  equations  (34)  and  (35),  X represents  the  axial  distance  from  the 
point  at  which  heat  or  mass  transfer  begins.  The  dimensionless  boundary- 
layer  thickness  6h+  is  the  same  in  both  equations  when  the  Prandtl 
number  equals  the  Schmidt  number  inasmuch  as  C+  = t+  for  a given  value 
of  y+  and  Prandtl  or  Schmidt  number. 

Values  of  local  Nusselt  and  Reynolds  numbers  can  be  obtained  from 
equations  (23),  (24),  and  (25),  as  for  fully  developed  flow,  with  the 
exception  that  the  expression  for  tt+  in  equation  (26)  is  replaced  by 


A similar  expression  can  be  obtained  for  Cb4-.  The  integral  in  the 
numerator  is  broken  into  two  parts,  because  t+  is  constant  and  equal 
to  tg+  outside  the  thermal  boundary  layer.  Inside  the  thermal  bound- 
ary layer  the  relation  between  t+  and  y+  is  obtained  from  figure  2. 
Values  of  u+  are  obtained  from  figure  1 for  y*  from  0 to  rQ+ 

inasmuch  as  a fully  developed  velocity  distribution  is  assumed.  The 
relation  between  Nusselt  number  and  X/D  for  various  values  of  Reynolds 
number  is  obtained  by  assuming  values  of  the  parameters  rQ+  and  5^+ 
and  by  calculating  the  various  quantities  from  equations  (23),  (25), 

(35),  and  (36). 

The  variation  of  local  Nusselt  number  for  heat  or  mass  transfer 
divided  by  the  fully  developed  Nusselt  number  with  x/D  and  Reynolds 
number  at  Prandtl  or  Schmidt  numbers  between  1 and  3000  is  given  in 
figure  6.  At  the  higher  Reynolds  numbers  the  values  of  Nu/Nud  for  a 

given  x/D  decrease  with  increasing  Prandtl  number;  that  is,  the  effect 
of  x/D  becomes  small  at  large  values  of  Prandtl  number.  At  low  Rey- 
nolds numbers  the  variation  is  more  complex:  Values  of  Nu/Nu^  first 

decrease  and  then  increase  slightly  as  Prandtl  number  increases.  In 
either  case  it  is  evident  that  in  the  entrance  region  the  fully  developed 
Nusselt  numbers  should  be  multiplied  by  a factor  which  is  a function 
of  X/D,  Reynolds  number,  and  Prandtl  number  rather  than  of  X/D  alone. 
That  is,  a simple  factor  such  as  (x/d)_c,  which  is  often  used  for  Prandtl 
numbers  of  approximately  1,  is  inadequate  for  high  Prandtl  numbers. 

The  effect  of  Reynolds  number  on  Nusselt  number  in  the  entrance  region 
increases;  that  is,  the  separation  of  the  curves  for  various  Reynolds 
numbers  increases  with  Prandtl  number.  The  same  conclusions  apply  to 
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the  average  Nusselt  numbers  plotted  in  figure  7 except  that  the  changes 
with  X/D  near  the  entrance  are  more  gradual  and  the  separation  of  the 
curves  with  Reynolds  number  is  greater  than  for  the  local  values. 

The  average  Nusselt  numbers  were  calculated  from 


Nuav  - 


x/d 

75 


Nu 


This  equation  is  consistent  with  the  definition 


r - tb)ctx 

(fc0  ■ ^av-  “ Y 


which  is  the  usual  way  of  defining  the  average  difference  between  wall 
and  bulk  temperatures  for  uniform  heat  flux. 

It  is  of  interest  to  note  that  the  Nusselt  numbers  for  turbulent 
flow  in  figures  6 and  7,  in  general,  display  trends  opposite  to  those 
for  Nusselt  numbers  for  laminar  flow  with  increasing  values  of  Prandtl 
number.  In  the  case  of  laminar  flow,  the  value  of  Nu/Nud  at  a given 

X/D  and  Reynolds  number  increases  considerably  with  Prandtl  number 
(ref.  26,  fig.  3),  because  the  heat  diffuses  through  the  fluid  more 
slowly  at  the  higher  Prandtl  numbers  (the  thermal  diffusivity  Is  lower) 
so  that  the  thermal  boundary  layer  is  thinner  and  Nu/Nu^  for  a given 
X/D  near  the  entrance  and  a given  Reynolds  number  is  consequently 
higher  than  for  the  lower  Prandtl  numbers.  The  same  phenomenon  also 
tends  to  increase  the  effect  of  X/D  for  turbulent  heat  transfer.  In 
the  case  of  turbulent  heat  transfer,  however,  the  shape  of  the  tem- 
perature profile  in  the  thermal  boundary  layer  changes  considerably  with 
Prandtl  number  and  becomes  very  flat  at  high  Prandtl  numbers  (fig.  2). 
This  means  that  the  temperature  profiles  for  fully  developed  flow  do 
not  differ  greatly  from  those  near  the  entrance  (both  are  flat)  although 
the  boundary-layer  thickness  for  the  two  cases  differs  considerably. 

The  Nusselt  numbers  in  the  entrance  region  for  turbulent  heat  transfer 
at  high  Prandtl  numbers  therefore  tend  to  quickly  approach  the  fully 
developed  values. 
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FULDf  DEVELOPED  HEAT  TRANSFER  AND  FRICTION  WITH  VARIABLE  PROPERTIES 

Turbulent  heat  transfer  to  liquids  vlth  variable  viscosity.  - In 
the  case  of  heat  transfer  to  liquids,  the  variation  of  the  viscosity 
with  temperature  is  considerably  greater  than  the  variation  of  the 
other  properties.  A good  approximation  to  the  actual  heat  transfer  in 
liquids  can  therefore  be  obtained  by  considering  only  the  viscosity  to 
vary  with  temperature.  Under  that  assumption  and  assumptions  1 and  2 
from  the  section  Analysis  for  Constant  Fluid  Properties,  equations  (4) 
and  (5)  become 


1= 

\ Mq  MO/PO  / dy+ 

and 

1 _ / J_  e \ 

\ Pr0  Mq/p 0 / dy+ 


(37) 


(38) 


For  variable  viscosity,  the  expression  for  e close  to  the  wall  (eq. 
(14))  can  be  written  in  dimensionless  form  as 


e 

Mo/pQ 


2 + + 
n u y 


nVy*  \ 
m/mo 


(39) 


For  liquids,  including  water,  oil,  ethylene  glycol,  and  sodium  hydroxide, 
u/n0  can  usually  be  represented  by  (t/tp)  , if  the  liquid  is  noi  too  near 
the  freezing  point.  The  exponent  d varies  from  -1  to  -4,  and  the  temper- 
atures are  measured  in  °F.  This  differs  from  the  case  for  gases  where  the 
temperatures  were  measured  in  °R  (ref.  l).  All  the  results  up  to  this 
point  are  independent  of  whether  the  temperatures  are  in  °R  or  °F.  From 
the  definitions  of  t+  and  3, 


u/mo  = (1  - pt+)d 


(40) 


Substituting  equations  (39)  and  (40)  into  equations  (37)  and  (38)  and 
writing  the  result  in  integral  form  yield 


(41) 
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and 


Equations  (4l)  and  (42)  can  be  solved  simultaneously  by  iteration,  that 
is,  assumed  values  for  u+,  y+,  and  t+  are  substituted  into  the  right 
side  of  the  equations  and  new  values  of  u+  and  t+  are  calculated  by 
numerical  integration.  These  new  values  are  then  substituted  into  the 
right  side  of  the  equations  and  the  process  is  repeated  until  the  values 

of  u+  and  t+  corresponding  to  each  value  of  y+  do  not  change 
appreciably.  Equations  (41 ) and  (42)  give  the  relations  between  u+, 
t+,  and  y+  for  various  values  of  the  heat-transfer  parameter  3 and 
of  Prg  for  the  region  close  to  the  wall. 

In  the  region  at  a distance  from  the  wall,  the  terms  in  the  equa- 
tions containing  variable  viscosity  are  neglected  so  that  u+  and  t+ 
are  given  by  equations  (21)  and  (22). 

Generalized  temperature  and  velocity  distributions  for  a Prandtl 
number  at  the  wall  of  10  and  a d of  -4  are  shown  for  various  values 
of  the  heat -transfer  parameter  3 in  figure  8.  The  value  y+  at  the 
intersection  of  the  curves  for  flow  close  to  and  at  a distance  from  the 

wall  is  taken  at  y^+  = 26,  as  in  the  case  of  constant  viscosity.  The 
effect  of  various  assumptions  for  the  variation  of  y^+  is  investigated 
in  reference  13  (fig.  13),  in  which  it  is  concluded  that  this  assumption 
should  give  accurate  results.  Positive  values  of  3 correspond  to  heat 
addition  to  the  liquid,  negative  values  to  heat  extraction.  The  values 
of  t+  at  a given  value  of  y*  increase  with  Increasing  3,  whereas 
the  values  of  u+  decrease.  These  opposite  trends  can  be  explained 
by  examining  the  terms  in  the  denominators  of  equations  (4l)  and  (42). 
The  first  term  in  the  denominator  of  equation  (4l)  (molecular)  causes 
u+  to  decrease  with  increasing  values  of  3 (d  < 0),  whereas  the 
second  term  (turbulent)  causes  u+  to  increase.  The  effect  of  the 
second  term  is  somewhat  smaller  so  that  the  net  effect  is  a decrease,  as 
shown  in  figure  8(b).  In  equation  (42),  however,  the  first  term  is 
independent  of  p and  therefore  the  second  term  causes  an  increase  in 
t+  with  increasing  values  of  p.  The  opposite  trends  for  heat  transfer 
and  friction  in  most  of  the  succeeding  curves  can  be  attributed  to  the 
same  cause. 
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The  effect  of  the  factor  In  parenthesis  in  equation  (42),  which  is 
equal  to  F in  equation  (13),  increases  with  increase  in  Prandtl  number 
because  of  the  steep  temperature  gradients  at  high  Prandtl  numbers  in 
the  region  very  close  to  the  wall  where  F differs  considerably  from  1. 
The  separation  of  the  t+  curves  for  various  values  of  3 (not  shown) 
therefore  increases  as  Prandtl  number  increases,  the  separation  for  a 
Prandtl  number  of  1 being  very  small. 

For  obtaining  the  relation  between  Nusselt  number,  Reynolds  number, 
Prandtl  number,  and  friction  factor  for  variable  viscosity,  equations 
(23),  (25),  and  (26)  apply  for  variable  properties  as  well  as  for  con- 
stant properties  if  the  viscosities  in  the  Reynolds  and  Prandtl  numbers 
are  evaluated  at  the  wall  temperature . The  friction  factor  f can  be 
M calculated  from 

& 

f = 2/(ub+)2  (43) 

Nusselt  numbers  and  friction  factors  are  plotted  against  Reynolds 
numbers  in  figures  9 and  10,  respectively,  with  the  viscosity  in  the 
Reynolds  and  Prandtl  numbers  evaluated  at  the  wall  temperature.  The 
curves  for  cooling  and  heating  are  for  values  of  db/po  on  order  of 

0.5  and  2,  respectively.  As  in  the  case  of  temperature  and  velocity 
distributions,  the  trends  with  increasing  values  of  3 for  the  Nusselt 
numbers  and  friction  factors  are  opposite.  Also,  as  in  the  case  of  the 
temperature  distributions,  the  separation  of  the  Nusselt  number  curves 
with  3 increases  as  the  Prandtl  number  increases,  the  separation  being 
very  small  for  a Prandtl  number  of  1.  Thus,  for  a Prandtl  number  of  1 
the  reference  temperature  for  evaluating  the  viscosity  in  order  to 
eliminate  the  effects  of  variable  viscosity  is  close  to  the  wall  tem- 
perature; the  departure  of  the  reference  temperature  from  the  wall  tem- 
perature increases  with  Prandtl  number.  This  is  to  be  expected  because, 
as  mentioned  previously  in  this  section,  the  separation  of  the  t+ 
against  y+  curves  with  3 increases  with  Prandtl  number.  In  those 
curves  all  the  properties  are  evaluated  at  the  wall  temperature . 

The  reference  temperature  for  heat  transfer  in  liquid  with  variable 
viscosity  for  a Prandtl  number  of  1 differs  from  that  obtained  for 
gases  in  reference  1,  in  which  all  the  properties  except  the  specific 
heat  were  considered  variable.  In  reference  1,  the  reference  temperature 
for  gases  was  found  to  be  close  to  the  average  of  the  wall  and  bulk 
temperatures  rather  than  close  to  the  wall  temperature.  In  assigning  a 
reference  temperature,  it  is  therefore  important  to  consider  what 
properties  are  variable. 

Values  of  x for  calculating  the  reference  temperature  tx,  where 
tx  * x(tQ  - tb)  + t^j,  are  given,  for  heat  transfer  and  friction,  in 


r 

i 


18 


NACA  TN  3145 


figure  11.  The  curves  for  3=0  in  figures  3 and  10  can  he  used  for 
variable  viscosity  if  the  viscosities  in  the  Reynolds  and  Prandtl  num- 
bers are  evaluated  at  the  reference  temperatures  given  in  figure  11 
(Rex  and  Prx).  The  values  of  x were  computed  for  values  of  d 
(h/mq  * (tAo)d)  an(d  “4  and  ?or  values  of  ^Aq  of  about  0.5 

and  2.  The  value  of  d had  little  effect  on  the  curves,  but  different 
curves  are  obtained  for  heating  and  cooling  of  the  liquid.  In  the  case 
of  heat  transfer,  the  reference  temperature  does  not  depart  greatly  from 
that  in  the  widely  used  Colburn  equation  (ref.  27),  wherein  the  viscosity 
is  evaluated  at  tn  ^ except  at  the  lower  Prandtl  numbers.  The  values 

of  x for  friction  are  lower  than  those  for  heat  transfer.  Deviations 
from  the  curves  in  figure  11  might  occur  for  very  high  viscosity  ratios 
or  for  cases  in  which  the  viscosity  variation  with  temperature  could 
not  be  represented  by  a simple  power  function. 

Turbulent  heat  transfer  to  gases  with  variable  properties  using 
present  method  of  analysis.  - If  the  present  method  of  analysis  is  to 
be  considered  more  general  than  the  analysis  in  reference  1,  it  must 
be  applicable  to  gases  with  variable  properties  as  well  as  to  liquids. 
When  assumptions  (l)  and  (2)  in  the  section  Analysis  for  Constant  Fluid 
Properties  are  used,  equations  (4)  and  (5)  become,  for  gases, 


where  Cp  is  assumed  constant  because  its  variation  with  temperature  is 
slight  compared  with  the  variations  of  viscosity,  thermal  conductivity, 
and  density.  As  in  reference  1,  it  is  assumed  that  k/kg  = = (tAp)d 

and  p/Pq  * tgA,  where  d s 0.68  and  the  temperatures  are  measured  in 
degrees  Rankine.  By  substituting  the  expression  for  e close  to  the 
wall  (eq.  (14))  and  tAo  = 1 - 3t+,  equations  (44)  and  (45)  can  be 
written  in  integral  form  for  the  region  close  to  the  wall  as 


dy+ 


n2uV 


8t+ 


n2u+y+  \ 

(1  - pt+)d+1  j 


(46) 
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Equations  (46)  and  (47)  can  be  solved  simultaneously  by  iteration  as 
were  equations  (41)  and.  (42). 

In  the  region  at  a distance  from  the  wall,  the  molecular  shear- 
stress  and  heat-transfer  terms  in  equations  (44)  and  (45)  are  neglected 
and  the  expression  for  e given  in  equation  (15)  is  used.  The  inte- 
gration is  carried  out  in  reference  1,  in  which  it  is  found  that 


where  y^+  = 26. 

Generalized  temperature  and  velocity  distributions  as  calculated  by 
the  present  analysis  for  a Prandtl  number  of  0.73  are  plotted  in  figure 
12.  The  distributions  calculated  in  reference  1 are  also  included  for 
comparison.  The  agreement  between  the  two  methods  of  analysis  is 
satisfactory,  so  that  the  Nusselt  numbers  and  the  reference  temperatures 
obtained  from  the  present  analysis  should  also  agree  with  those  given 
in  reference  1. 

Laminar  heat  transfer  to  liquids  with  variable  viscosity  and 
uniform  heat  flux.  - Laminar  heat  transfer  to  liquid  metals  with  variable 
viscosity  is  investigated  in  reference  28.  The  results  of  that  inves- 
tigation should  be  applicable  also  to  liquids  with  high  Prandtl  numbers 
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if  the  value  of  d is  on  the  order  of  -1.6.  The  temperatures  in  the 
report  are  to  be  measured  in  °F  (u/m-q  = (tAn)d)*  It  Is  found  (ref. 
28)  that  the  results  for  fully  developed  laminar  flow  with  constant 
heat  flux  can  be  represented  closely  by 

>••(8)0°'“ 


SUMMARY  OF  RESULTS 

The  following  results  were  obtained  from  the  analytical  investi- 
gation of  heat  and  mass  transfer  in  smooth  tubes  at  high  Prandtl  and 
Schmidt  numbers: 

1.  By  modifying  the  expression  for  eddy  diffusivity  from  a previous 
analysis  to  account  for  the  effect  of  kinematic  viscosity  in  the  region 
close  to  the  wall,  good  agreement  was  obtained  between  predicted  and 
experimental  results  for  heat  and  mass  transfer  at  Prandtl  and  Schmidt 
numbers  between  0.5  and  3000.  A simplified  equation  was  obtained  for 
very  high  Prandtl  or  Schmidt  numbers. 

2.  The  analysis  indicated  that,  except  at  low  Reynolds  numbers, 
the  entrance  effect  (local  Nusselt  number  divided  by  fully  developed 
Nusselt  number)  for  heat  or  mass  transfer  decreases  as  Prandtl  or 
Schmidt  number  increases. 

3.  The  analysis  indicated  that  the  effects  of  variable  viscosity 
on  turbulent  heat  transfer  and  friction  in  liquids  can  be  nearly 
eliminated  in  ordinary  cases  by  evaluating  the  viscosities  in  the  Rey- 
nolds and  Prandtl  numbers  at  reference  temperatures  which  are  functions 
of  the  Prandtl  number.  For  the  laminar  case  the  results  for  liquid 
metals  with  variable  viscosity  given  in  a previous  analysis  should  be 
applicable  to  liquids  at  high  Prandtl  numbers. 

4.  When  the  present  method  of  analysis  was  applied  to  gases  with 
variable  properties,  essentially  the  same  results  were  obtained  as  are 
reported  from  analysis  and  experiments  in  a previous  report. 
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APPENDIX  - SYMBOLS 

The  following  symbols  are  used  in  this  report: 

concentration  of  diffusing  substance,  (lb) (sec2)/ft4 

bulk  concentration  of  diffusing  substance,  (lb) (sec2)/ft4 

concentration  of  diffusing  substance  at  y ■ 5, 

(lb) (sec2)/ft4 

concentration  of  diffusing  substance  at  wall, 

(lb) (sec2)/ft4 

exponent 

specific  heat  of  fluid  at  constant  pressure,  Btu/(lb)(°F) 

specific  heat  of  fluid  at  constant  pressure  at  wall, 
Btu/(lb)(°F) 

inside  diameter  of  tube,  ft 

exponent,  value  of  which  depends  on  variation  of  viscosity 
of  fluid  with  temperature 

function  of  e'/(p/p) 

acceleration  due  to  gravity,  32.2  ft/sec2 

local  heat-transfer  coefficient,  ^/(^q  - t-^), 

Btu/ (sec) (sq  ft)(°F) 

local  mass-transfer  coefficient,  n^/(CQ  - C^),  ft/sec 

average  heat-transfer  coefficient,  %/ - t^)av 

average  mass-transfer  coefficient,  n^/(tQ  - )av- 

thermal  conductivity  of  fluid,  Btu/ (sec) (sq  ft)(°F/ft) 

thermal  conductivity  of  fluid  evaluated  at  tQ, 

Btu/ (sec) (sq  ft)(°F/ft) 

rate  of  mass  transfer  toward  tube  center  per  unit  area, 
(lb) (sec)/cu  ft 
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rate  of  mass  transfer  toward  tube  center  per  unit  area 
at  wall,  (lb)(sec)/cu  ft 


constant 

rate  of  heat  transfer  toward  tube  center  per  unit  area, 
Btu/(sec)(sq  ft) 

rate  of  heat  transfer  at  inside  wall  toward  tube  center 
per  unit  area,  Btu/(sec)(sq  ft) 

inside  tube  radius,  ft 

temperature,  °F  for  liquids  or  °R  for  gases 

bulk  static  temperature  of  fluid  at  cross  section  of  tube, 
°F 

reference  temperature  for  local  Reynolds  and  Prandtl 
numbers,  x(tg  - t-^)  + t^,  °p 

temperature  of  fluid  outside  thermal  boundary  layer,  °F 

wall  temperature,  °F  for  liquids  or  °R  for  gases 

average  difference  between  wall  and  bulk  temperature,  °F 

time-average  velocity  parallel  to  axis  of  tube,  ft/sec 

bulk  velocity  at  cross  section  of  tube,  ft/sec 

distance  from  point  at  which  heat  or  mass  transfer 
begins,  ft 

number  used  in  evaluating  arbitrary  temperature  in  tube 
*x 

distance  from  wall,  ft 

thermal  or  diffusion  boundary -layer  thickness,  ft 

coefficient  of  eddy  diffusivity  for  momentum,  sq  ft/sec 

expression  for  eddy  diffusivity  which  neglects  effect  of 
M-/ p , sq  ft/sec 

coefficient  of  eddy  diffusivity  for  heat  or  mass,  sq  ft/sec 


3159 


NACA  TN  3145 


x Karman  constant 

X molecular  diffusivity,  sq  ft/sec 

Xq  molecular  diffusivity  at  wall,  sq  ft/sec 

|i  absolute  viscosity  of  fluid,  (lb)(sec)/sq  ft 

Mb  absolute  viscosity  of  fluid  evaluated  at  t , 

(lb)(sec)/sq  ft  “ 

|i  absolute  viscosity  of  fluid  evaluated  at  t , 

(lb)(sec)/sq  ft  X 

Mn  absolute  viscosity  of  fluid  evaluated  at  tn, 

(lb)(sec)/sq  ft 

p mass  density  of  fluid,  (lb) (sec^)/ft4 

PQ  mass  density  of  fluid  evaluated  at  tQ,  (lb) (sec2)/ft 

•t  shear  stress  in  fluid,  lb/sq  ft 

Tq  shear  stress  in  fluid  at  wall,  lb/sq  ft 


Dimensionless 


groups : 

concentration  parameter. 


bulk  concentration  parameter, 


(Cp  ~ cb )\f To/pQ 

“o 


f 


Nu 

Nu' 

Nu. 


av 


Nu 


av 


value  of  C+  at  &h+,  (CQ  ~ C5^  /Tp/p0 

“o 

value  of  C+  at  y^+ 
friction  factor,  2Tq/pu^^ 

Nusselt  number  for  heat  transfer,  hD/k 
Nusselt  number  for  mass  transfer,  h'D/x 
average  Nusselt  number  for  heat  transfer,  h&vD/k 
average  Nusselt  number  for  mass  transfer,  hav'D/X 


t 
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fully  developed  Nusselt  number  for  heat  transfer 
fully  developed  Nusselt  number  for  mass  transfer 
Peclet  number,  pu^Dcpg/k 
Prandtl  number,  cpg|i/k 

Prandtl  number  with  properties  evaluated  at  tg 
Reynolds  number,  pu^D/p 

Reynolds  number  with  viscosity  evaluated  at  tx,  pu^D/i^ 
Reynolds  number  with  viscosity  evaluated  at  tg,  pu^D/|ig 

w To/pn 

tube  radius  parameter,  Vg/p  r0 
Schmidt  number,  u/(pk) 

Schmidt  number  at  wall,  Mq/po\) 

Stanton  number  for  heat  transfer,  h/pgubcp 
Stanton  number  for  mass  transfer,  h'/u^ 

temperature  parameter,  -^-2 — ■ ^ P,^0-  = ^ 0 

^ Vp0  ^ 

bulk- temperature  parameter,  i 1 “ T~  ) 

P \ *0  / 


»(-£) 


V 

value  of  t+  at  y + 

y 

u+ 

velocity  parameter,  u/  / 

V 

bulk -velocity  parameter 

V 

value  of  u+  at  y^+ 

y To/po 

V^o/po~ 


VP0 


wall  distance  parameter. 
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value  of  y+  at  intersection  of  curves  for  flow  close  to 
wall  and  at  a distance  from  wall 

ratio  of  eddy  diffusivities,  e^/e 

heat-transfer  parameter,  q^ *\J tq/ p^j  (CpS^o^ 

dimensionless  thermal  or  diffusion  boundary- layer  thick- 
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Figure  2.  - Generalized  temperature  or  concentration  distributions  for  various  Prandtl  or  Schmidt  numbers. 
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Figure  4.  - Comparison  of  analytical  and  experimental  results  for  fully  developed  heat  and  mass  transfer 


Concluded.  Comparison  of  analytical  and  experimental  results  for  fully  developed  heat  and  mass 


Figure  5.  - Comparison  of  various  analyses.  Reynolds  number,  10,000. 
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.gore  7.  - Variation  of  average  Nusselt  number  over  fully  developed  Nusselt  number  for  heat  or  mass  transfer  with  X/D 
and  Reynolds  number.  Uniform  wall  heat  flux,  uniform  initial  temperature  distribution,  and  fully  developed  velocity 
distribution. 
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Concluded.  Generalized  distribution  for  liquid  with  variable  viscosity.  Fully  developed  flow;  Prandtl 
iAi0  = (t/t0)~4.  (Vertical  line  is  dividing  line  between  eqs.) 
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Figure  12.  - Comparison  of  generalized  distributions  for  gases  with  variable  properties.  Prandtl  number,  0.73.  (Vertical 
line  is  dividing  line  between  eqs.) 
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Figure  IE.  - Concluded.  Comparison  of  generalized  distributions  for  gases  with  variable  properties.  Prandtl  number 
0.73.  (Vertical  line  is  dividing  line  between  eqs.) 


